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Abstract 

Disorder along a DNA strand due to non uniformity associated with the counter ion type and 
location, and in rise and twist are investigated using density functional theory. We then model the 
conductance through a poly(G) DNA strand by including the influence of disorder. We show that 
the conductance drops by a few orders of magnitude between typical lengths of 10 and 100 nm. 
Such a decrease occurs with on-site potential disorder that is larger than 100 meV. 


PACS numbers: 73.63.-b, 85.65. +h, 87.15.-v, 72.80. -r 
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I. INTRODUCTION 


During the last few years, the possibility to use DNA in nanotechnology and as a nano- 
device has been increasingly investigated. This interest stems from the ability of controlled 


growth of nucleotide sequences and DNAs surprising conducting properties. Transport mea- 
surements 1-7 have led to a broad range of conductances, opening a huge controversy on this 
issue. The value of DNA conductance would be very sensitive to various factors some of 
which are the nature of the metal contacts, defects, conformational changes and the location 
of the Fermi energy. In this paper, we model the role of counter ions along the back bone of 
dry DNA in causing conductance reduction with increase in strand length. The distance be- 
tween base pairs in B-DNA is approximately 3.4 A and the hybridization energy (U) between 
the HOMO levels of adjacent guanines of a poly(G) sequence is small (around 87 meV). 3 ’ 8 ’ 9 
The corresponding HOMO band width (4 U) is a narrow 350 me’V . Small changes in the 


environment of +Vin 
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pairs or irregularities along the length of the poly(G) 
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bound to exist in dry DNA. For example, upon drying, it is unlikely that the counter ions 
of the phosphate group will all be identical and periodically arranged. It is also conceivable 
that a counter ion could be missing from a phosphate group. Such non idealities will lead to 
change in the energy eigen values (on-site potential) of base pairs along the strand. If this 
perturbation is small compared to U, then transport through poly(G) will be unaffected. 
The main idea of this paper is that the eigenvalues of base pairs are modified by values com- 
parable to U because of irregularities even for poly(G). Then electronic delocalization arising 
from the n stacking of the base pairs is broken over long distances, and the conductance 
through a poly(G) molecule will significantly diminish with increase in length. 

We investigate the influence of various structural irregularities such as change in rise 
and twist, and ion mismatches on the eigenvalues of a dimer of guanines. The ab initio 
calculations are performed using Gaussian 98 10 in the density functional approximation 
with the B3YLP 11,12 exchange-correlation functional. The basis set used here corresponds 
to 6-21G. More diffuse or complete basis sets like the 6-3lG(d) do not affect significantly 
our results. Each base pair is assumed to have only a single counter ion as discussed below . 
We note that work on dry DNA that clarifies the number of counter ions around each base 
pair would be useful and could quantitatively change our results. We have compared the 
eigenvalues corresponding to the HOMO, LUMO and the average difference between all 
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the eigenvalues of the system (Table 1). Change of rise and twist (by 6°) do not lead to 
very large variations of the HOMO and LUMO eigenvalues. Variation in the counter ion 
characteristics affect the HOMO and LUMO eigenvalues more significantly. The HOMO 
eigen value changes by 100 meV and 50 meV in the cases where both H 3 0 + are replaced by 
Na + and K+. Variations of the distance between the counter ion and the phosphate group 
also lead to variations of the order of 50 meV. A larger difference, nearly 500 meV and 1000 
meV, in the HOMO and LUMO values is found between hydrogen (bonded to one oxygen 
of the phosphate) and H 3 0 + and Na + counter ions. We should note that the occurrence of 
a hydrogen bonded to the phosphate is rather unlikely. Over all, the change in eigenvalues 
in table I makes our suggestion that irregularities along a poly(G) causes a reduction m 
conductance with length plausible. 

The conductance is calculated using the Landauer-Buttiker approach. 13 The Hamiltonian 
for an extended system of base pairs is contracted by starting from the Hamiltonian of 
dimers, following the nearest neighbor approximation scheme of references 14 and 15. The 
Hamiltonian of the dimers are calculated using Gaussian 98 10 and the transmission proba- 
bility is calculated using the scheme in reference 16. A key ingredient to the calculation is 
the strength of the disorder potential. The disorder is introduced via a uniform change in 
the on-site potential of all atoms in a base pair by 5e h where * corresponds to the base pair 
labeled i. Further, it is expected that for two adjacent nucleotides, there is some degree 
of correlation between the amplitude of disorder potential at neighboring base pairs. A 
correlation in the variation of the on-site potential of the base pairs is introduced using the 
following procedure. The disorder potential of the first nucleotide in the sequence is chosen 
randomly in the interval ± | e M |. The disorder of the second nucleotide is chosen randomly 
subject to the constraint that de 2 cannot differ from by more than *f. This procedure is 
continued for the other base pairs. A typical potential profile is shown in Fig. 1. To justify 
the use of the smoothing procedure, we have compared the eigenvalues of dimers with two 
Na + , two H 3 0 + and a Na + and H 3 0 + each. We find that the change in the HOMO level 
of the third case is in between the first two cases. Finally, in order to focus on the effect 
of disorder, the open boundaries correspond to disoreder-free semi infinite nanotube leads 

instead of realistic metal contacts. 
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For a DNA molecule trapped between two electrodes, the HOMO / LUMO are the most 
important in determining the conductance. Here, we investigate the conductance when the 
Fermi energy is close to the HOMO channel. 2 The system considered is a poly(G) without 
the backbone. This should be a good approximation because for regular B DNA the HOMO 
is due to the bases and the backbone only plays a role in the structural stability. The 
transmission and conductance (in inset) for a disorder-free poly(G) molecule is presented m 
Fig. 2. The HOMO channel is located at around -5.5 eV and has a width of approximately 
350 meV. At zero bias, the maximum possible conductance observable with this molecule 
is around 2G 0 (where G 0 c- i^fkn)- The experimentally observed conductances have been 
orders of magnitude smaller. Various factors such as DNA-contact coupling, conformational 
changes of the DNA lying on the metal and on the substrate will play a role in determining 
the conductance. Also, mismatch in the work function of the the metal and DNA wall pla> 
a role in determining the conductance. For example, in a long strand of DNA, the band 
bending between DNA regions lying on the metal contact and over the substrate could 
be larger than the HOMO band width. Neglecting the above mentioned issues (especially 
the last item which is important for systems with a small energy band width), we focus 
on the role of disorder in affecting the linear response conductance. Considering that the 
HOMO channel width is of the order of 350 meV, it is expected that the transmission 
will decrease drastically for devices longer than the localization length. 17 This is what we 
basically observe. The conductance (for a Fermi level of -5.55 eV which corresponds to the 
maximum value of the conductance obtained without disorder) as a function of the length of 
the device for different value of the disorder strength e M (is shown in Fig. 3). The exponential 
decrease over the length scale shown is not observed for small disorders (e M - 50 meV) and 
the conductance remains of the order of unity till lengths up to 300 nm (not presented). 
However, for value of e M larger than 100 meV, the conductance tends to drop very rapidly 
with length over the length scale of 100 nm considerered, due to the large mismatch of the 
on-site potential between the nucleotides. The plot shown for e M = 150 meV is an example 
of this. Furthermore, considering that small disorders (smaller than 100 meV) do not lead to 
significant drop of the conductance, over a length scale of lOOnm, it seems \ery unlikely that 
twist or rise mismatches (of the strength shown in Table 1) can explain the phenomenon of 
very large resistances observed in references 2 and 6. Mismatch in the counter ions could lead 
to a large decrease in conductance with length but such an effect can only be verified if both 
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short and long DNA strands are suspended between the metal contacts. An experimental 
setup such as in reference 18 may be able to verify this prediction. 

In conclusion, we have shown that disorder induced by variation in the specifics (type and 
distance) of the counter ions along a DNA strand is likely to hinder the delocalization of the tt 
orbitals even in a poly(G) molecule. We find that while random changes in the twist and rise 
affect the eigenvalues of the guanines, a non crystalline counter ion arrangement, can lead 
to significant changes in the eigenvalues of base pairs when compared to the hybridization 
energy between guanines in poly(G). We calculated the coherent conductance through a 
poly(G) and find that the conductance of strands over a length scale of 100 nanometers 
is significantly affected when counter ions cause a change in the eigenvalues of guanines 
by more than about 100 meV. Modeling of the coupling of electrons to vibrational degrees 
of freedom, which will lead to realistic values of the residual conductance with disorder is 

beyond the scope of this work. 
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A enoMO (meV) 

Aslumo (meV) 

Ae (meV) 

Na + vs. HsO" 5 " h 

93.4 

136.4 

125.9 

K 4- vs. H 3 0 + 

69.6 

89.8 

106.7 

Na + , Na + vs. Na+, H 3 0 + 

49.1 

65.3 

59.7 

H 3 0+, H 3 0 + vs. Na + , H 3 0+ 

44.3 

71.1 

68.1 

K + vs. Na + 

23.8 

46.6 

58.6 

H- vs. H 3 0 + 

552.8 

1054.4 

1161.1 

H- vs. Na + 

646.2 

1190.8 

1285.5 

Na + distance (Ad = +0.1 A) 

27.1 

51.0 

48.3 

Na + distance (Ad = —0.2 A) 

63.1 

119.4 

114.6 

twist (Afi = ±2°) 

2.0 

5. 

32.0 

twist (Afi = —6°) 

0.9 

19.5 

115.6 

twist (Af2 = +6°) 

16.9 

10.4 

73.1 

rise (Ar = —0.1 A) . 

17.4 

12.9 

47.7 

rise (Ar — +0.1 A) 

14.6 

11.1 

43.9 

rise (Ar = —0.2 A) 

37.9 

28.6 

98.2 

rise (Ar = +0.2 A) j 

26.6 

20.8 

81.5 


TABLE I: Eigenvalue differences (for the HOMO and LUMO eigenvalues and also for the average 
differences between all the eigenvalues) between various dimmers. Three kinds of system varia- 
tions are presented: Mismatch of the counter ions (first 7 rows), mismatch of the distance of the 
counterions with respect to the phosphate group (next 2 rows) and mismatch of the twist and rise 

of the dimmer (next 7 rows). 
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. 1: Amplitude of the disorders as function of the length of the device. The solid line corresponds 
maximum amplitude of the disorder of 50 meV and the dashed line to 150 meV 
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FIG. 2: Conductance versus fermi energy for an infinite poiy(G) DNA molecule in a B conformation 
without disorders. 
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FIG. 3: Conductance (at 
as function of its length, 
the maximum amplitude 


T = 298 


K and a Fernai .eve! of -5.55 eV) of a poly(G) B DNA mdecu.e 

and triangles correspond to cases where 


The lines with diamonds, stars 


of the disorder is 50 meV , 100 


meV and 150 meV respectively. 
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